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Abstract. We consider the dynamics of one or more self gravitating shells of matter in a 
centrally symmetric gravitational field in the Painleve family of gauges. We give the reduced 
hamiltonian for two intersecting shells, both massless and massive. Such a formulation is applied 
to the computation of the semiclassical action of two intersecting shells. The relation of the 
imaginary part of the space-part of the action to the computation of the Bogoliubov coefficients 
is revisited. 



1. Introduction 

In this paper we shall study the dynamics of one or more self-gravitating spherical shells of matter 
subject to a centrally symmetric gravitational field and the application of the ensuing formalism 
to the semiclassical treatment of Hawking radiation. Such field of research was started by the 
papers by Kraus and Wilczek [TJ [2]. The main results we shall report here are the extension 
of the treatment to one massive shell of matter and also the extension to two or more massive 
shells of matter which in the time development can also intersect [3]. The main appeal of the 
approach of [TJ [2] is that energy conservation is taken exactly into account and thus the back 
reaction effects can be computed. 

We shall keep the formalism for more that one shell as close as possible to the original 
formalism of [I] which can be summarized in words that we shall adopt a Painleve gauge. 

The extraction of the reduced hamiltonian is often classified as a very complicated procedure. 
We shall show that by introducing a "generating function" the derivation can be drastically 
simplified and also by the same token it can be extended to deal not only with massive shells 
but also with a finite number of massive shells. The treatment of massless shells is just a 
particular case. 

In Sec|2] we shall give the derivation of the reduced action in the case of one shell while 
proving the independence of the canonical momentum within the Painleve class of gauges. In 
Sec|3] we shall discuss the equations of motion. In SecH] we work out the analytic properties of 
the conjugate momentum p c which appears in the reduced action and in SecG2 we shall briefly 
describe the extension of the treatment to more than one shell. An application of the results 
is given in SecE] giving a derivation within such a formalism of the Dray-'t Hooft- Redmount 
exchange relations. An important integrability relation is reported in Seed which allows to 
compute the imaginary part of the space integral of the canonical momenta for two shells. In 
SecJS]we revisit the role of the imaginary part of the canonical momentum in determining the 



Bogoliubov coefficients and thus the most important features of Hawking radiation. In SecEJwe 
give some concluding remarks. 



2. The reduced action 

As usual we write the metric for a spherically symmetric configuration in the ADM form 

ds 2 = -N 2 dt 2 + L 2 (dr + N r dt) 2 + R 2 dQ 2 . (1) 

were following [TJ [21 HI Ej we shall choose the functions N, N r , L, R as continuous functions of the 
coordinates. We shall work on a finite region of space time tt) x (ro, r m ) . On the two initial 
and final surfaces we give the intrinsic metric by specifying R(r,ti) and L(r, tj) and similarly 
R(r,tf) and L(r,tf). 

The complete action in hamiltonian form, boundary terms included is [TJ ® [7] 

rtt rr m 

S = S sheU + dt dr(ir L L + tt r R - NH t - N r H r ) + 

Jti Jro 

rt f . , . NRR' 



dt (-N r iT L L 
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where 

SsheU = dtpf. (3) 

Jti 

r is the shell position and p its canonical conjugate momentum. 7i r and Tit are the constraints 

n r = 7r R R' -7r' L L-pd(r-f), (4) 

JU ~ RR " + + L7Tl RR ' L ' ™ L + J r?L-* + m* Sir f) (5) 

Action ([2]) is immediately generalized to a finite number of shells. The shell action as given 
by eqs. (|3!4l5p refers to a dust shell even though generalizations to more complicated equations 
of state have been considered El [9]. The boundary terms are those given in the paper by 
Hawking and Hunter [6] and will play a very important role in the following. The function F [3] 




2-M / R ' ( R '\ 2 , , 

-R- +RRlos [T-nT) " 1 + irl (6) 

has the remarkable property of generating the conjugate momenta as solutions of the constraints 
as follows 

nL = lL = dL (7) 
_ 6F_ _ dF_ _ d_dF_ 
nR ~ 6R~ dR~ frdR 1 ^ ' 

where A4 is a mass which is constant in r except at the shell position f and that we shall denote 
by H for r above all the shell positions and by M for r below all the shell positions. We shall 
adopt a Painleve gauge defined by L = 1 everywhere. With regard to the remaining freedom 
in the choice of the gauge we shall choose R = r except for a deformation region near the shell 
positions. Such deformation is unavoidable because the constraints impose a discontinuity in 
the derivative R'(r) at r = f as we shall see in the following. In the Painleve gauges F becomes 

F = RW(R, R', M) + RR'(C(R, R',M) - B{R, M)) (9) 



where 



and 



W(R,R',M) = Jr' 2 "i + ^r; C{R,R',M) =log(R' -W(R,R',M)) (10) 



B(R,M) = J^ + log\ 1 



' 2M 
~R 



(11) 



where we exploited the freedom of adding to ([6]) a total derivative, thus gaining for F the useful 
property of vanishing wherever R' = 1. For L = 1 eqs. ()7l8p become 



vr L = RdR' 2 - 1 + — = W; 7tr = — . 

With regard to R(r, t) one can choose several gauges, within the Painleve family. 

R R 
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Figure 1. Outer gauge 



Figure 2. Inner gauge 




Typical are the "outer gauge" characterized by R(r, t) = r for r > f(t) and shown in FigfTJ 
and the inner gauge characterized by R{r,t) = r for r < fit) shown in FigJ21 but there are also 



more general choices as shown in Figj3j We stress that such deformation g is not related to the 
thickness of the shell which is always zero. 

Using the solutions (|7|8p of the constraints one obtains in the outer gauge the following 
reduced [3] action where only the shell position f appears as degree of freedom 



J (Per- M(t) J —dr + (-ATttz + NRR%™ dt. (13) 



rt f ( i r (i) dF 

'ro 

The p c is easily computed 



Pc = f{AC - AB) (14) 

where AC = C{f+e)—C{r—e) and similarly for AB. In deriving eq. (|14p we used the consequences 
of the constraints (|4|5p 

V 



AR' = — — ; where V = yp 2 + m 2 ; Aitl = —p. (15) 
R v 



In the outer gauge we find 



Pc = V2M~r -Vm-f-r log ( + ^ ^ 



with p given implicitly by 



m 2 I2H 



H - M = V + — -p\ /— ; V = \ p 2 + m 2 . (17) 
2r y r v 

This is the result obtained by Friedman, Louko, Winters-Hilt [4J through a limit procedure in 
which the support of the deformation function goes to zero but actually it is independent of the 
deformation. For working out the dynamics of one shell the limit procedure is all right, but in 
dealing with two or more shell it is important to have a deformation on a finite range because 
otherwise a limit procedure would give a result dependent on the order in which the two limits, 
deformation and f\ — > r~2, are taken. 

Similarly one can compute using the general formula the reduced canonical momentum of 
the system in the inner gauge 

pi = V2Mr - V2Hr-rlog ( - ^ 2Hr ^ \ ( 18 \ 

\f - V + p- V2Mf J 

and p, again determined by the discontinuity equation (fT5|) . is given now by the implicit equation 



m 2 „ I2M 



H-M = V p\ ; V = Jp 2 + m 2 . (19) 

2r y r v 

The two reduced canonical momenta p c , p % c appear to be completely different but they can be 
proven to be the same function of f. One can consider also the more general gauges depicted in 
Fig. 3 and it can be proven [10] that p c is always the same. However in inner gauge a term H(t) 
appears in the reduced action while in the more general gauges both a term M and a term H 
appear in the action. 

The boundary term given in eq.© is equivalent to 



-HN(r m )+MN(r ). 



(20) 



3. Equations of motion 

In deriving the equation of motion from action (|13|) one can consider M(t) = M, the interior 
mass, as a datum of the problem and vary H to obtain 



*im " N{rm) = (21) 

which using the expression of p c and the relation between N and iV r imposed by the gravitational 
equations 0] can be written as 



f = ^iV(f) - N r (r) =11-^] N(f). (22) 

Alternatively one can consider H, the total energy as a datum of the problem and vary M(t). 
The calculation is far more complicated due to the presence of M in the action (|13p . but using 
the equations for the gravitational field one reaches [3] the same equation of motion (|22p . In 
addition a consequence of the gravitational equations of motion is the constancy in time of M(t) 
and H(t). 

4. Analytic properties of p c 

We saw that in the outer gauge p c is given by (|16ll7p . The solution of eq. (fT7|) for p is 



(23) 



where 



H — M m 2 

A = — is- (24) 

If we want p to describe an outgoing shell we must choose the plus sign in front of the square 
root. Moreover the shell reaches r = +oo if and only if H — M > m as expected. 

The logarithm in p c , eq.(|16|). has branch points at zero and infinity and thus we must 
investigate for which values or r such values are reached. At f = 2H, p has a simple pole 
with positive residue. Thus the numerator in the argument of the logarithm in p c goes to zero 
and below 2H it becomes 

f-V-p- V2Hf (25) 



where here V is the absolute value of the square root. Expression (|25p is negative irrespective 
of the sign of p and stays so for f < 2H because p is no longer singular. As a consequence p c 
below 2H acquires the imaginary part mf. Below r = 2M the denominator of the argument 
of the logarithm in eq. (|16p becomes negative so that the argument of the logarithm reverts to 
positive values. Thus the so called classically forbidden region where p c becomes complex is 
2M < r < 2H independent of m and of the deformation g and the integral of the imaginary 
part of p c for any deformation g and for any mass m of the shell is 

r f 2H A S 
Im / p c dr = vr / rdr = 2tt(H 2 - M 2 ) = (26) 

J J2M 2 

which is the original result derived in [TJ [11] in the zero mass case. Parikh and Wilczek [11] gave 
to 

exp(-2Im J p c dr) (27) 
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Figure 4. Two shell dynamics 



the interpretation of the tunneling probability for the emission of a quantum of energy uj. 
Criticism and alternative proposals for the emission probability like 

exp(— Im <j> p c dr) (28) 

exp ^— 2Im( J p c dr + temporal contribution)^ (29) 

followed [T21 (13 [Til 1151 HZl US H3l I2JH |2Sl EDI EB ■ Here instead we shall discuss the role of 
eq. p^j) in the framework of mode analysis. This will be done in Sec®. 

5. Two shell reduced action 

Now we have three characteristic masses, M, H and an intermediate mass Mq, see Fig. 4, which 
can change only if the two shells cross. Working as before with M = const considered as a 
datum of the problem we reach the reduced action [3] 

dT ■ dT 

hPd + f 2 Pc2 + H(R{h) - h)^V + M (R(h) - + 

d r 2 „ , ■ n dF 



where 



+ dtJ ro m dr + (-N^ L + NRR% 7 (30) 



V2_ 



T = log^ 7 f T ; V = R(AC - AB)\ fl ; (31) 



p cl = R'(h+e)V; p c2 =p° c2 + —(R^ 1 )-f 1 )V (32) 

dr 2 

and p^ C 2 is given by eq. ()16p with M replaced by Mq. The novelty is that now even in the outer 
gauge the time derivative of H intervenes in addition to Mq and p c \ and p c2 depend both on 
fx and f 2 . We can vary r%, f 2l H and Mq independently obtaining the correct equations of 
motion [3]. As expected one finds that the exterior shell moves irrespective of the dynamics 
which develops at lower values of r until a crossing occurs. 




Figure 5. The exchange diagram 

6. The exchange relations 

In case of crossing of the two shells from the equations of Sec.Q we can obtain relations between 
M, Mq, H,r e and Mq being r e the the shell position at the crossing and Mq the intermediate 
mass after the crossing. During the crossing the masses of the shells can change, provided they 
satisfy a relation analogous to the energy-momentum conservation is special relativity 

Pi+P2= p'i + P 2 (33) 

Vl + V 2 = V[ + V^ V n = ^pl + ml; = V 1 p'l + m*. 

This is an outcome of the constraints. It is not possible to predict the final masses of the two 
shells as it depends on the details of the interaction which has to be specified. A relatively 
simple case is when the masses are unchanged during the crossing (transparent crossing) [3] and 
the simplest case is the crossing of two massless shells which remain massless, thus re-obtaining 
the well known Dray-'t Hooft-Redmount relations [211 [22] which we report here below 

Hr e + Mr e - 2HM = M f e - 2MqM' q + M' Q f e (34) 

being f e the crossing radius. 

7. Integrability of the form p c idr\ + p c2 df 2 and imaginary part of the space- 
component of the action 

The space-part of the on-shell action for two massive shells of matter is given from eq. (|30p by 

/ (Pd fi +p c2 r 2 ) dt (35) 
Jti 

as on the equations of motion H{t) = M{t) = 0. It is possible to prove a theorem analogous to 
the one found in the books of Whittaker and Arnold [231 [HI i- e - that in presence of a constant of 
motion, in addition to the hamiltonian, the form p c \df\ +p C 2dr 2 is closed even though the proof 
is somewhat different [3]. The intermediate mass Mq plays the role of the additional constant 
of motion. The above result allows to deform the integration path as to bring f\ immediately 
below or above f 2 - This is allowed by the absence of discontinuities in our scheme. But, in 
words, two coaleshed shell have the same properties of a single shell with energy H — M. The 
final result is that 

Im / dt(p cl ri +p c2 f c2 ) = Im / ' (p cl dr x +Pc2 df 2 ) = 2tt(H 2 - M 2 ). (36) 

The result (|36p holds also for in the case when at the crossing the two shells can change their 
mass. For details we refer to [3]. 



8. Mode analysis 

The original way to extract information on the spectrum of the radiation is mode analysis 
[26], [U [27] . In the massless case 

J r p c df' = f(f,M)-f(f,H) (37) 
can be computed exactly and its expansion to first order in co = H — M is 

/f 
p c df = 4Mw log(r — 2M) + regular terms. (38) 

Given the semiclassical mode 

( j ) {r,t)=e i r p ^'- iujt (39) 

we can perform the analysis of the mode regular at the horizon in terms of the above mode. 
This analysis can be performed either by scalar product [25] i.e. space integration where one has 
to keep into account that the background metric is Painleve or by time Fourier analysis. In this 
way one obtains the well known formulas for the Bogoliubov coefficients. The regular outgoing 
modes near the horizon are given by 

^f,t) = e ^-2M )e -OT_ (40) 

Computing the scalar product, and taking into account that the background metric is the 
Painleve metric, we have 

- i J (rd p - 0d p iP*)g pO eore<f,V=g dfdOd^ (41) 

with 

g rt = N r = ^M ^ (42) 

and the integration region is outside the horizon. Keeping only the most singular terms at the 
horizon eq. (j4"Tj) reduces, with r = exp(— t/4M), to 

g— ikxr ^AiMu) log x— iujt (43) 
X ' 

We can compute such scalar products at t = obtaining the standard Hawking integral 

/ e -ifcx e 4*Af w logx ^_ = e 2^M nA-^ M TUiuM) = const at, (44) 
./() x ...... 

which gives the dominant contribution for large k. The coefficient (3^ is obtained by changing 
in (|44|) to into —u. Alternatively we can extract the Bogoliubov coefficients by performing a 
time Fourier transform i.e. 

/+oo r+oo r j T 

g— ikxTgiiMu log x— iu)t ^ / g— ikxr ^AiMu) log(xr) (45) 
-oo JO T 

which is, as it should be, a result independent of r (i.e. x) and reproduces eq.(|44p. The derivation 
is valid to first order in uj which is the realm of validity of the external field approximation. For 



finite oj it is problematic to perform a space integration on modes, because we have not a well 
defined background metric. 

Thus Kraus and Wilczek |lj followed the time Fourier analysis method. In order to do so one 
has to construct the non-perturbative modes regular at the horizon. This is not a completely 
trivial task. The regular modes are given by 

e iS = e ikr(0)+i /; (0) P cdf-i{H-M)t ^ 

where S is the on shell action computed with the following boundary conditions: 1) At time t 
the shell position is r, outside the horizon. 2) At time t = the conjugate momentum is k, i.e. 
S(0,f) = kr. Due to these boundary conditions r(0), H and as a consequence p c depend on 
k, t, r even if along each trajectory H is always a constant of motion. Using the saddle point 
approximation [TJ [27] one obtains that the absolute value of the Bogoliubov coefficient is 
given by 

■ iQ )Pc (r>,H,M)dr> = ^ImJ^ Pc (r',H,M)dr^ 

computed for H = M + u. What is important here is that only the space part of the action 
appears, as the time part (H — M)t cancels with tot at the saddle point. r(0) is given by the 
condition 

k = r m o g ^E-^E (48) 

which for H = M + uj > M is solved by 

2M < 2H < f (0) (49) 

and thus there is no imaginary contribution to the integral from the gap 2M, 2H as we discussed 
in Sect.dU). The Bogoliubov coefficient (3^ is obtained by changing in (f4"T|) uj in —uj (always 
uj > 0). Then we have 

r(0) < 2H < 2M (50) 
The value of the saddle point t, now complex, is given by 




(51) 

VH0)~V2H 

and the r.h.s. term of eq. (|47p becomes 

e -4TrMu)(l-u)/2M) 

where the last passage is due to eq. ()26p for which we gave a general proof within the family 
of Painleve gauges even in the massive case. An explicit derivation of the result (|52p has been 
given in [10] by working out the late time expansion of the time development of the action. The 
main point is that in the mode analysis what comes in is not the total action of a model particle 
crossing the horizon, but only the "space part" of it. The analysis and the results depend on 
the interpretation of the expression 

exp (iS) (53) 

as modes of the field dressed by the gravitational interaction with non vanishing quanta of 
energy. The result eq. (|36p proven for the emission of two shells [3] was interpreted [28] as the 
absence of correlations among quanta in the emitted radiation. It would be of interest to give 
a similar "mode interpretation" of the result (|36p derived for the emission of two shells which 
during the time evolution can also interact PJ. To this end one should compute the two shell 
modes which are regular at the horizon, and perform a time Fourier analysis of them. This has 
not yet been accomplished. 



9. Concluding remarks 

In this paper we gave a general treatment of the dynamics of one or more self- gravitating 
spherical shell of matter in a spherical gravitational field. We extended the treatment of pQ to 
more than one shell [3] proving on the way the universality of the reduced conjugate momentum 
within the family of Painleve gauges [ lfjj . 

This allows to give a general derivation of the exchange relations and also, exploiting an 
integrability result, to extend the result on the imaginary part of the space part of the action 
to more than one shell. Instead of following the tunneling picture here we have revisited the 
treatment which follows by interpreting the exponential of the action, properly subtracted, as 
the dressed modes of the Hawking radiation. The main point here is that according such a mode 
interpretation only the space part of the action and in particular its imaginary part comes in 
determining the Bogoliubov coefficient and thus in determining the Hawking spectrum corrected 
for the back reaction effects. 
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